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Preface 


The  structural  effects  of  the  outer  scale  of  turbulence  are  significant  in  certain 
propagation  problems  in  the  Earth’s  lower  atmosphere.  In  particular,  image 
wander  (angle  of  arrival  fluctuations  of  distant  objects)  and  shape  distortion 
of  objects  due  to  nonisoplanatism  are  critically  dependent  on  the  shape  of 
the  refractive  index  power  spectrum  in  its  energy-containing-eddy  region.  To 
correctly  model  these  effects  requires  an  accurate  model  of  the  outer  scale  of 
turbulence,  both  as  a  function  of  height  and  as  regards  how  it  impacts  the 
refractive  index  spectrum. 

In  the  past,  the  outer  scale  was  defined  rather  vaguely,  using  order  of  magnitude 
arguments.  Or,  at  best,  outer  scale  was  defined  based  of  the  mathematical  form 
used  in  describing  the  refractive  index  spectrum  itself.  Also,  the  outer  scale  was 
often  equated  with  the  so-called  mixing  length  used  in  numerical  weather  model 
closure  methods  for  solving  the  Navier-Stokes  equations. 

But  the  mixing  length  is  merely  a  convenient  order-of-magnitude  length 
scale.  Atmospheric  optics  has  tended  to  adopt  this  value  in  lieu  of  better 
parameterizations.  This  report  attempts  to  fill  this  data  void  by  dealing  with 
several  aspects  of  the  outer  scale  problem.  To  accurately  determine  outer 
scale  effects  with  regard  to  propagation  problems  involving  optical  turbulence 
(turbulent  image  wander,  image  shape  distortion,  etc.)  better  models  must  be 
sought  than  those  based  on  reasonable  guesses. 

First,  I  propose  a  definition  for  the  outer  scale  which  is  independent  of  the 
particular  mathematical  form  used  in  describing  the  refractive  index  spectrum. 

Second,  an  analysis  of  temperature  spectrum  data  from  the  1968  Kansas 
experiment  is  discussed  that  leads  to  a  new  form  for  the  refractive  index 
spectrum.  For  completeness,  I  also  propose  a  combined  spectrum  called  the 
composite  spectrum  that  includes  both  inner  and  outer  scale  influences. 

Third,  based  on  the  analysis  of  the  spectrum  itself,  a  relationship  is  determined 
between  the  peak  frequency  data  from  the  Kansas  experiment  and  the  outer  scale 
of  turbulence.  This  relationship  describes  a  functional  form  for  the  turbulence 
outer  scale  as  a  function  of  both  height  above  the  surface  and  a  surface  layer 
stability  parameter,  the  Monin-Obukhov  length. 
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Executive  Summary 

Synopsis 

This  report  considers  the  nature  of  the  outer  scale  of  turbulence,  its  relationship 
to  the  refractive  index  power  spectrum,  and  the  form  of  this  spectrum  based  on 
data  reported  from  the  1968  Kansas  experiment  [Kaimal  et  al.,  1972]. 

The  outer  scale  of  turbulence  is  a  length  scale  used  to  demarcate  the  beginning 
of  the  so-called  inertial  subrange.  For  two  points  separated  by  distances  greater 
than  the  so-called  inner  scale  length  (a  few  millimeters)  and  less  than  the  outer 
scale  length  (on  the  order  of  meters)  fluctuations  in  refractive  index  follow  a 
“2/3  law”.  Difficulties  arise  for  separation  distances  outside  this  range.  These 
deviations  are  the  subject  of  spectral  modifications  that  account  for  both  inner 
and  outer  scale  variations. 

Under  many  circumstances  the  influence  of  the  outer  scale  is  not  significant. 
For  example,  optical  turbulence  effects  on  beam  propagation  generally  involve 
scintillation  and  beam  spread  influences  that  are  relatively  impervious  to  outer 
scale  effects  due  to  the  limited  width  of  the  beam.  Imaging  of  point  sources 
through  turbulence  is  also  little  affected  by  outer  scale  when  considering  only 
turbulent  blurring  and  scintillation  effects  since  the  optics  entrance  aperture  is 
smaller  than  the  outer  scale. 

However,  for  object  jitter  and  image  distortion  effects,  outer  scale  influences 
are  likely  dominant.  The  outer  scale  length  has  been  shown  to  be  directly 
related  to  the  spatial  frequency  that  has  the  greatest  influence  on  angle  of  arrival 
variations  [Tofsted,  1992].  Yet  outer  scale  is  not  well  defined  or  quantified. 
To  correct  deficiencies  regarding  outer  scale  effects  characterization,  this  report 
focuses  on  three  issues:  (1)  Defining  outer  scale  independent  of  the  particular 
mathematical  form  being  adopted  to  describe  the  refractive  index  spectrum. 
(2)  Developing  a  more  accurate  refractive  index  spectral  form  for  low  spatial 
frequency  characterization.  (3)  Deriving  an  equation  for  evaluating  outer  scale 
based  on  meteorological  data  measurements. 

Results 

To  define  the  outer  scale  of  turbulence  requires  a  knowledge  of  the  form  the 
refractive  index  spectrum  assumes  when  inner  and  outer  scale  effects  are  not 
considered.  This  spectrum  is  called  the  Kolmogorov  spectrum, 

$«*(*)=  0.033  C^k-11/3,  (1) 
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where  k  is  a  spatial  frequency  variable  given  in  units  of  radians  per  meter  and 
Cn  the  refractive  index  structure  parameter  with  units  of  refractive  index 
squared  (unitless)  times  meters  to  the  -2/3  power. 

The  Kolmogorov  spectrum  is  the  most  frequently  encountered  spectrum  in  the 
optical  turbulence  literature  because  it  has  the  simplest  form.  This  property 
makes  the  Kolmogorov  spectrum  a  useful  benchmark,  as  follows:  Assume  we 
identify  some  refractive  index  spectral  form  which  accounts  for  outer  scale  effects. 
Call  it  $„(«).  For  k  values  within  the  so-called  inertial  subrange  we  will  have 
$«(«)  ~  ®uk(k).  For  k  below  this  range,  one  enters  the  energy-containing- 
eddy  region  in  which  <&n  <  with  the  ratio  $>n/®nK  approaching  zero  as  k 
approaches  zero.  Within  this  framework,  the  outer  scale  can  be  defined  as, 


Though  arbitrary,  heretofore  no  quantitative  definition  has  been  established. 
Using  this  definition,  the  Kansas  1968  experimental  data  reported  by  Kaimal  et 
al.  [1972]  were  analyzed.  While  the  goal  of  the  analysis  was  to  evaluate  the 
outer  scale  length,  it  was  necessary  to  first  produce  a  new  model  of  the  refractive 
index  spectrum  for  low  spatial  frequencies,  k.  This  was  necessary  because  the 
Kansas  data  were  plotted  as  one-dimensional  (1-D)  spectra  {Fn(n i)),  typical 
of  sensor  data  collected  through  meteorological  tower  measurements,  where  the 
wind  advects  spatially  fluctuating  air  past  the  tower  location.  On  the  other  hand, 
atmospheric  turbulence  studies  are  based  on  three-dimensional  (3-D)  spectra. 
The  only  way  to  compare  the  two  is  to  propose  a  3-D  spectrum  and  transform 
it  to  its  1-D  representation,  as, 


OO 

?n(«l)  =47 T  J  \  +  Kr 


(3) 


Comparison  of  spectra  led  to  the  development  of  a  relationship  between  outer 
scale  and  the  reported  peak  frequency  of  a  universal  function  of  scalar  property 
fluctuations  in  turbulent  flow.  This  peak  frequency  is  given  in  terms  of  z/L , 
the  dimensionless  ratio  of  height  above  the  surface  (z)  to  the  Monin-Obukhov 
length  ( L ). 


Conclusion 

Together,  the  new  definition  of  outer  scale  developed,  the  new  spectral  form 
to  describe  the  effects  of  the  outer  scale  on  the  refractive  index  spectrum,  and 
the  new  relationship  describing  outer  scale  as  a  function  of  height  and  stability 
extend  the  accuracy  of  modeling  low  spatial  frequency  effects  dramatically. 


1.  Introduction 


The  atmospheric  optics  community  has  generally  relied  on  a  few  select  forms  for 
representing  refractive  index  spectra,  as  illustrated  by  the  following  examples. 
Nelson  et  al.  [2000]  use  a  simple  Kolmogorov  spectrum  [Kolmogorov,  1962] 
that  exhibits  n~n/3  behavior  at  all  spatial  frequencies,  where  k  is  a  radially 
symmetric  spatial  frequency  variable  given  in  units  of  radians  per  meter.  (This 
spectrum  was  described  in  the  Executive  Summary.)  Frehlich  [2000]  augments 
the  standard  Kolmogorov  spectrum  to  account  for  inner  scale  effects  through 
the  addition  of  a  Hill  [1978]  Spectrum  term  which  accounts  for  the  peculiar 
bump  phenomenon  which  occurs  in  the  vicinity  of  the  frequency  equivalent  of 
the  inner  scale  length.  These  studies,  and  others  [Frehlich,  2000;  Davis  and 
Walters,  1994;  Martin  and  Flatte,  1988;  etc.],  focus  on  beam  wave  propagation, 
the  goal  of  which  is  usually  to  limit,  as  much  as  possible,  the  horizontal  spread 
of  such  beams.  The  narrow  cones  used  in  these  studies  tend  to  not  experience 
significant  outer  scale  effects.  And  the  phase  screens  used  in  propagating  the 
beams  are  usually  centered  on  the  beam  path. 

Perhaps  more  sensitive  to  outer  scale  effects  are  studies  such  as  Ellerbroek  [1997] 
and  Young  et  al.  [1998].  The  former  uses  a  von  Karman  spectrum  referenced 
to  Tatarskii  [1971].  The  latter  uses  a  modified  von  Karman  spectrum  which 
includes  an  inner  scale  roll  off  region  generally  attributed  to  Tatarskii  [1971]. 
Similarly,  Yan  et  al.  [2000]  used  a  von  Karman  spectrum  in  their  adaptive  optics 
study. 

Together,  the  Kolmogorov  spectrum,  in  combination  with  von  Karman, 
Tatarskii,  and/or  Hill  spectral  modifiers,  generally  span  the  gamut  of 
the  available  choices  of  spatial  frequency  spectra  for  characterizing  optical 
turbulence.  Significant  theoretical  and  experimental  work  has  been 
accomplished  in  the  high  spatial  frequency  regime,  with  various  representations 
of  this  bump  effect  available  [Churnside,  1990;  Andrews,  1992;  Frehlich,  1992]. 

All  spectra  exhibit  similar  k-11/3  dependence  for  k  within  the  inertial  subrange, 
and  the  behavior  of  the  Andrews  and  Frehlich  models  are  very  similar.  Of  more 
concern  are  differences  between  proposed  versions  of  the  von  Karman  spectrum. 
While  all  models  of  this  spectrum  are  given  as, 


$n(«) 


0.033  Cl 
(k2  +  k  I) 11/6 


(1) 


there  is  disagreement  over  the  definition  of  k©  as  well  as  problems  with  the 
behavior  of  this  functional  form  at  low  spatial  frequencies. 
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For  example,  Ishimaru  [1978],  Ellerbroek  [1997]  and  Young  et  al.  [1998]  define 
Kq  as  1/La,  where  La  is  referred  to  as  the  outer  scale.  But  Kopeika  [1998] 
and  Clifford  [1978]  define  ft©  as  2tt / L0.  yet  also  refer  to  L0  as  the  outer  scale. 
Obviously  two  divergent  forms  of  the  same  equation  cannot  lead  to  identical 
definitions  for  the  outer  scale. 

Ishimaru  indicated  La  should  be  on  the  order  of  the  height  above  the  surface  in 
appendix  C  of  his  book.  Kopeika  set  La  to  kv  z,  where  kv  is  von  Karman’s 
constant,  equal  to  approximately  0.385  over  rough  surfaces  [Frenzen  and 
Vogel,  1995],  and  2  is  the  height  above  the  surface.  Combining  the  uncertainty 
about  the  definition  of  La  with  the  uncertainty  about  the  definition  of  ft©,  we 
have  over  an  order  of  magnitude  disagreement  between  Ishimaru  and  Kopeika 
regarding  the  positioning  of  the  outer  scale  cutoff. 

The  von  Karman  spectrum  also  has  a  second  problem.  At  the  low  frequency 
limit,  as  ft  — >  0,  the  spectrum  should  also  approach  a  zero  value.  The 
von  Karman  spectrum  approaches  the  constant  value  0.033  C2  ft©11^3  instead. 
The  zero  limit  requirement  can  be  shown  through  the  relationship  between  the 
refractive  index  spectrum  $n  and  the  spatial  autocorrelation  function,  Tn(r) 
[Goodman,  1985],  given  by, 

T„  (r)  =  (n1  ( x )  n'  (x  +  f)) .  (2) 

Here,  the  angle  brackets  represent  an  expectation  operation,  n'  represents 
refractive  index  perturbations  at  positions  x  and  x+r  about  an  average  refractive 
index  value. 

$n  is  related  to  Tn  through  a  Fourier  relationship  as  will  be  discussed  further 
in  section  2.3.  However,  because  T„  is  an  autocorrelation  function  of.  n',  <hn 
must  be  proportional  to  the  magnitude  squared  of  the  Fourier  transform  of  n' 
[cf.,  Gaskill,  1978].  But  n'  is  a  zero  mean  quantity  by  definition  and  hence 
(Parseval’s  theorem)  nn s  Fourier  transform  must  equal  zero  at  the  origin.  By 
extension,  <hn  must  also  equal  zero  at  the  origin.  Neither  the  Kolmogorov  nor 
the  von  Karman  spectra  exhibit  this  property. 

Any  study  purporting  to  consider  optical  turbulence  effects  on  image 
propagation  should  begin  by  considering  the  outer  scale  influences  on  the 
problem  because  the  outer  scale  represents  the  size  scale  with  the  greatest 
fluctuation  energy  in  the  spectrum.  The  available  energy  normally  can  be 
measured  as  a  function  of  <kn(ft)  times  4itk2  dn,  the  differential  volume  associated 
with  fluctuations  at  a  particular  spatial  frequency  scale.  Since  is  decreasing 
as  ft  /  across  the  inertial  subrange,  the  combined  effect  will  be  the  maximum 
fluctuations  that  occur  somewhere  near  the  outer  scale  size. 


2.  Analysis 


To  address  the  deficiencies  in  the  existing  theory  discussed  in  the  introduction, 
several  aspects  of  optical  turbulence  related  to  low  spatial  frequency  phenomena 
must  be  considered.  A  new  form  for  the  refractive  index  spectrum  is  proposed 
to  replace  the  currently  dominant  von  Karman  spectrum.  Spectral  analysis 
arguments  are  used  to  indicate  the  limiting  behavior  of  the  spectrum  at  low 
spatial  frequency.  Analysis  of  actual  data  provides  further  insight  into  the 
detailed  behavior  of  the  spectrum  in  the  vicinity  of  the  outer  scale  frequency, 
1/L0.  However,  because  of  ambiguities  in  definition  of  the  outer  scale  itself, 
it  is  necessary  to  clarify  the  meaning  of  this  quantity  first  and  then  use  this 
more  precise  definition  in  relating  the  outer  scale  to  measureable  meteorological 
quantities  associated  with  the  atmospheric  surface  layer. 

In  section  2.1,  a  new  definition  of  outer  scale  is  proposed.  This  new  definition 
is  linked  to  a  property  of  the  outer-scale-influenced  refractive  index  spectrum 
itself,  rather  than  to  a  specific  form  used  in  describing  that  spectrum.  In 
section  2.2,  the  use  of  the  mixing  length  is  briefly  discussed  that  has  led  to  some 
confusion  in  the  past.  In  section  2.3,  an  analysis  is  performed  to  link  the  research 
findings  for  the  Kansas  experiment  to  the  three-dimensional  (3-D)  refractive 
index  spectrum.  The  results  of  this  analysis  permit  both  postulating  a  new 
outer  scale  influenced  spectrum  and  establishing  a  relationship  between  outer 
scale  and  the  Monin-Obukhov  length,  an  atmospheric  surface  layer  stability 
parameter.  The  implications  of  this  functional  behavior  on  the  refractive  index 
itself  will  then  be  considered,  along  with  a  discussion  in  section  2.4  of  a  composite 
spectrum  that  accounts  for  both  inner  and  outer  scale  effects. 

2.1  Outer  Scale  Defined 

The  first  step  in  describing  the  outer  scale  effects  should  be  to  uniquely  identify 
the  outer  scale  in  such  a  way  that  it  is  not  dependent  on  the  form  of  the 
equation  in  which  it  appears.  As  was  discussed  in  chapter  1,  the  appearance  of 
the  outer  scale  in  the  von  Karman  spectrum  results  in  an  ambiguity  based  on 
different  literature  sources.  The  same  value  of  outer  scale,  La,  has  significantly 
different  ramifications  when  defining  «©  =  1/L0  [Ishimaru,  1978;  Andrews 
and  Phillips,  1998]  as  opposed  to  the  use  of  kq  =  2n /L0  [Clifford,  1978; 
Kopeika,  1998]. 
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To  resolve  this  dilemma  I  propose  the  following:  Let  represent  the 

Kolmogorov  spectrum, 


$»*(«)  =  0.033  C*K-11/3,  (3) 

where  n  is  a  spatial  frequency  variable  with  units  of  radians  per  meter  and  C%  is 
the  refractive  index  structure  parameter  with  units  of  refractive  index  squared 
(unitless)  times  meters  to  the  -2/3  power.  The  refractive  index  spectrum  thus 
has  units  of  volume.  Let  $n(/«,  £0)  represent  any  spectrum  where  outer  scale 
effects  are  included,  and  where  the  parameter  £0  reflects  outer  scale  effects  but 
is  not  necessarily  equal  to  the  outer  scale.  To  quantitatively  define  the  value  of 
the  outer  scale  length,  L0,  we  need  a  new  variable,  k0  =  l/L0.  We  will  define 
the  outer  scale  as  that  length  which  satisfies  the  relation, 

(ftQ)  ■£(?)  _  1  ,  .v 

$nK  (k<z)  2 

Typically,  £„  appears  in  functional  forms  describing  <f>n  involving  only  the 
product  kC0.  Due  to  the  way  £„  appears  in  these  functional  forms,  outer  scale 
parameters  in  these  equations  can  always  be  written  as  linear  relations  relative 
to  the  actual  outer  scale:  £0  —  aL0. 

For  example,  in  the  case  of  the  von  Karman  spectrum,  when  k0  =  l/£0  is 
used,  we  must  have  a  equal  to  (26/11  —  l)  ^  ~  1.475.  Otherwise,  using 

the  definition,  kq  =  2t r/£0  a  is  written  as,  (26/11  -  l)-1/2  »  9.269.  In 
either  case  La  maintains  the  same  position  relative  to  the  knee  of  the  spectrum 
regardless  of  the  form  adopted  for  kq. 

2.2  Outer  Scale  versus  Mixing  Length 

One  hazard  to  avoid  in  the  available  literature  is  the  difference  between  outer 
scale  and  a  quantity  customarily  called  the  mixing  length  or  the  characteristic 
length.  While  various  mixing  length  definitions  are  available,  these  are 
usually  employed  in  providing  closure  mechanisms  for  solving  the  Navier  Stokes 
equations  [e.g.,  Lewellen,  1977;  Zhao,  1994;  Finardi  et  a/.,  1995].  In  this 
literature  Kopeika’s  kvz  value  often  appears.  For  example,  Lewellen  defines 
his  characteristic  length  as, 


A0 


1.68 


du/dz 


1.68 


kyZ 


4m(z/L)' 


(5) 


where  it*  is  the  friction  velocity,  du/dz  is  the  vertical  gradient  of  the  horizontal 
windspeed,  and  (j)m  is  a  diabatic  influence  function  of  momentum  which  is  a 
dimensionless  function  of  the  ratio  of  the  height  (z)  divided  by  the  Monin- 
Obukhov  length  ( L ).  I  previously  used  this  parameter  for  defining  outer  scale 
[Tofsted,  1992],  though  in  hindsight  it  appears  inappropriate.  Zhao  [1994], 


12 


and  Finardi  et  al.  [1995],  use  a  form  of  the  mixing  length  scale  developed  by 
Blackadar  [1962]: 


1  _  1  1 

Aqq  hyZ 


(6) 


where  Aqo  is  the  asymptotic  value  of  A0  at  large  height.  However,  in  no  instance 
is  this  result  linked  to  a  turbulence  spectrum.  Hence,  applications  requiring  an 
outer  scale  should  look  to  other  sources  to  determine  this  quantity. 

Such  a  source  seems  to  be  exclusively  available  through  data  illustrated  in 
or  derived  from  Kaimal  et  al.  [1972].  Nevertheless,  use  of  this  information  is 
complicated  by  the  mode  of  presentation  of  the  data.  Optical  turbulence  work 
usually  involves  integration  of  3-D  spectra,  but  the  Kansas  experimental  data 
presented  by  Kaimal  involves  one-dimensional  (1-D)  spectra.  Hence,  to  use  the 
Kaimal  data  a  turbulence  spectrum  must  be  posited,  and  then  this  spectrum 
must  be  converted  to  a  1-D  representation.  This  task  is  treated  in  the  following 
section.  This  reference  can  also  be  found  along  with  many  of  the  other  articles 
cited  in  this  report  in  the  SPIE  Milestone  series  compilation  -  Turbulence  in  a 
Refractive  Medium,  Andreas  [1990]. 


2.3  One-  and  Three-Dimensional  Spectra 

Within  the  optical  turbulence  research  community  the  assumption  of  3-D 
homogeneous  and  isotropic  turbulence  is  typically  made.  This  assumption  is 
useful,  though  in  most  cases  it  is  false.  Hinze  [1987]  points  out  that  homogeneous 
turbulence  is  only  possible  under  decaying  turbulence  conditions.  To  ensure 
that  the  turbulence  is  maintained  requires  nonhomogeneous  energy  sources. 
Nevertheless,  under  most  conditions  the  source  terms  generating  additional 
turbulence  are  small.  I  thus  adopt  a  locally  homogeneous  condition. 

Theory  also  suggests  that  turbulence  is  not  isotropic  at  k  <  1/La.  However, 
again,  while  use  of  the  isotropic  assumption  is  not  entirely  appropriate,  the 
largest  deviations  occur  only  at  the  lowest  frequencies,  and  these  lowest 
frequencies  have  the  least  effect  on  the  optical  properties  normally  studied. 
Therefore,  from  a  tractability  view  and  due  to  its  minimal  impact,  an  isotropic 
turbulence  assumption  will  be  adopted  throughout  this  discussion.  In  general, 
though,  we  do  not  need  to  adopt  a  homogeneous  turbulence  assumption  in 
order  to  study  properties  of  the  refractive  index  power  spectrum.  Further,  it  is 
generally  assumed  that  turbulence  has  a  vertically  inhomogeneous  structure. 
With  these  caveats  made,  I  proceed  with  the  general  development,  first 
considering  the  limitations  of  the  von  Karman  spectrum. 

As  was  mentioned  in  the  introduction,  the  refractive  index  spectrum  is  related 
to  the  autocorrelation  function,  Tn(r).  Or,  more  generally,  for  scalar  property 
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£,  the  relationship  between  its  autocorrelation  function  and  its  spectrum  may 
be  written  as  [Beland,  1993], 

3  00 

$£  («)  =  (2^)  JJJ  r€  (r)  exp  (* « •  f)  df,  (7) 

—  OO 


where  k  =  («i,  «2,  «3)>  r  =  (#,  2/,  2),  and  i  =  V— 1* 

The  Fourier  relation  between  and  permits  the  writing  of  the  inverse 
transform  relation, 


—  OO 


(«)  exp  (— i  k  •  f)  d,K, 


(8) 


When  the  assumption  of  isotropy  is  invoked  these  equations  may  be  reexpressed 
in  spherical  coordinates  and  integrated  over  their  spherical  angular  coordinates 
first,  leaving  functions  dependent  only  on  the  magnitudes  k  =  |k|  and  r  =  |f|. 

3  00 

*»<«)-(£)  /r«(r)^4^*,  (9) 

0 

00 

(r)  =  f  (k)  s^n(Atr)  47rK2  (10) 

J  ft  T 

0 


Most  atmospheric  optics  calculations  are  made  in  terms  of  weighted  integrals 
involving  <Fn.  But  the  refractive  index  is  a  function  of  temperature,  pressure, 
and  humidity  fluctuations  in  the  atmosphere,  with  temperature  fluctuations 
being  the  dominant  effect  for  visible  to  infrared  frequencies.  To  measure 
these  fluctuations  directly  usually  involves  tower  measurements  which  record 
the  temporally  fluctuating  properties  of  the  atmosphere  as  the  wind  advects 
positionally  varying  air  through  the  sensor  mast  location.  Taylor’s  [1938]  frozen 
turbulence  hypothesis  is  then  used  to  relate  the  temporal  fluctuations  observed 
to  spatial  fluctuations  through  extrapolation  based  on  the  mean  wind  speed. 

The  temporal  records  of  the  tower  data  can  be  analyzed,  but  these  data  only 
represent  1-D  spectra.  These  spectra  are  related  to  their  3-D  counterparts 
through  the  relation  [Panofsky  and  Dutton,  1984], 

OO 

Ft  («l)  =  2  J  $£  27 TKrdKr,  (11) 

0 

where  nr  is  a  polar  spatial  frequency  coordinate  in  the  (k2,  K3)  plane, 
corresponding  to  spatial  coordinates  normal  to  the  wind  direction.  Note  also 
that  Fjf  is  only  defined  for  positive  K\  values.  Due  to  function  symmetry 
about  zero,  results  at  negative  K\  values  are  folded  onto  the  positive  portion 
of  the  axis.  This  approach  leads  to  the  factor  of  2  in  front  of  the  integral  in 
eq.  (11). 


2.3.1  Atmospheric  Refractive  Index 

I  now  note  that  although  f  was  used  to  denote  an  arbitrary  scalar  property  of 
the  atmosphere,  a  so-called  conservative  passive  additive,  our  real  interest  is  in 
fluctuations  of  the  atmospheric  refractive  index.  Optical  turbulence  is  caused 
by  fluctuations  in  this  index  about  its  mean  value.  The  mean  refractive  index 
is  a  function  of  air  temperature,  pressure,  and  humidity.  Positional  changes  in 
these  properties  alter  the  refractive  index  of  the  atmosphere  (n). 

In  fact,  n  is  so  close  to  unity  in  the  atmosphere  that  a  secondary  variable,  N, 
the  refractivity,  is  normally  used  to  describe  n’s  variability: 

n  =  1  +  JV  x  10-6.  (12) 


The  functional  dependence  of  N  is  given  in  terms  of  its  dry  air  behavior  in 
combination  with  water  vapor  effects  [Hill  et  al.,  1980], 


N  —  Ndry  -|-  -/V vapor ,  (13) 

P  N 

Ndry  =  0.3789  x  [1  +  (5.337  -  0.0157  T)  Pt  x  10“6]  (14) 

where  Pt  is  the  air  pressure  in  torr,  and  T  is  the  air  temperature  in  Kelvin  (K). 
No  expresses  the  wavelength  dependence,  given  by, 


Nq  =  64.328  + 


29498.1 
(146  -  A"2)  + 


255.4 

(41 -A"2)’ 


(15) 


where  A  (n m)  is  the  radiation  wavelength  in  the  visible  through  far  infrared 
band. 

The  water  vapor  influence  on  the  refractive  index  is  given  by, 

Nvapor  =  -0.0059  x  (1  -  0.0109  A“2)  A,  (16) 

where  A  is  the  absolute  humidity  (g/m3). 

At  visible  through  far  infrared  wavelengths,  although  there  are  humidity  and 
pressure  dependencies,  temperature  variations  tend  to  clearly  dominate  the 
fluctuations  in  the  refractive  index  about  its  mean.  A  secondary  correlation 
between  temperature  and  humidity  fluctuations  has  been  characterized  by 
Wesely  [1976],  though  the  main  correlation  effect  can  be  accounted  for  through 
use  of  the  Bowen  ratio  [Kunkel  and  Walters,  1983].  These  temperature 
fluctuations  are  driven  largely  by  the  energy  flux  of  sensible  heat  at  the  earth’s 
surface  [Oke,  1978]. 

To  quantify  the  variations  in  refractive  index  due  to  temperature,  we  write  the 
differential  of  n  with  respect  to  temperature, 


(17) 


dn  —  dN  x  1(T6  »  -78.2  x  10"6  ^  dT, 
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where  pressure,  P,  has  been  converted  to  units  of  millibars  (mbar)  (1  torr  = 
1.33322  mbar).  The  important  aspect  to  be  emphasized  here  is  that  fluctuations 
in  refractive  index  are  directly  proportional  to  fluctuations  in  temperature.  We 
can  thus  equivalently  study  either  a  turbulent  field  of  temperature  fluctuations 
or  refractive  index  fluctuations,  and  we  will  be  studying  the  same  phenomenon. 

To  directly  study  these  fluctuations,  the  refractive  index  is  normally  rewritten 
as  an  average  quantity,  no,  and  a  fluctuation  component,  n' . 


n  =  n0  +  n'.  (18) 

There  are  no  optical  turbulence  effects  due  to  no,  while  the  n'  statistic  was 
previously  discussed  in  connection  with  the  autocorrelation  function,  Tn(r). 


2.3.2  Template  Three-Dimensional  Spectra 

Hence  temperature  and  refractive  index  fluctuations  are  strongly  coupled,  and 
once  the  structure  of  the  temperature  spectrum  is  understood,  the  refractive 
index  spectrum  immediately  follows,  having  the  same  form.  The  difficulty 
remains,  however,  in  determining  an  appropriate  3-D  spectrum  when  only  1- 
D  measurement  data  are  available. 


To  solve  this  problem  a  family  of  template  3-D  spectra  will  be  used.  These  can 
be  transformed  into  1-D  spectra  and  analyzed  with  respect  to  available  empirical 
curves  and  data  reported  in  the  literature.  This  family  of  functions  has  the  form, 


Tm  (k,  C)  = 


Cll/3  (/c£)2 m 
,  ,  2l  (H+em)/6 » 

1  +  (k£)2j 


(19) 


such  that  for  k  C,  YTO  (/c,  C)  — >  k-11/3.  Based  on  this  family,  the  standard 
von  Karman  spectrum  [Tatarskii,  1971]  of  turbulence  can  be  written,  using  the 
m  =  0  member  of  this  family,  as, 


$nV  («)  =  fin  C„Y o  («,  £o)  ■ 
where  (3n  is  an  exact  form  replacing  the  use  of  0.033, 


5  22/3r(5/6) 
36  7r3/2T(2/3) 


«  0.033. 


(20) 


(21) 


In  defense  of  Theodore  von  Karman,  I  would  note  that  it  is  arguable  that  he 
did  not  actually  propose  the  spectrum  honoring  his  name,  von  Karman  was 
certainly  aware  that  the  limiting  behavior  of  turbulent  spectra  at  low  frequencies 
approaches  a  zero  value  at  the  origin.  Hinze  [1987]  attributes  the  equivalent  of 
Ti  to  von  Karman  [1948],  while  he  attributes  the  use  of  To  to  Saffman  [1967]. 
All  members  of  this  family  behave  similarly  at  high  spatial  frequency  k,  where 
they  exhibit  k  11 7/3  dependence.  But  at  low  spatial  frequencies  To  approaches 
a  nonzero  constant. 
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2.3.3  The  Kaimal  Curve 

To  study  the  behavior  of  this  family  in  comparison  to  available  data,  the 
Tm  functions  were  first  transformed  to  1-D  spectra  using  eq.  (11).  These  1- 
D  functions  were  then  multiplied  by  frequency  variable  ki  so  they  could  be 
compared  to  a  ‘universal’  empirical  curve  suggested  by  Kaimal  et  al.  [1972]. 
This  equation  was  most  recently  expressed  by  Kaimal  and  Finnegan  [1994]  in 
the  form, 

-  l  W/°) 

O I  l  +  fl(///o)5/3’  1  ’ 

where  £  is  the  property  under  study  (temperature,  humidity,  etc.),  <r|  is  the 
variance  of  this  property,  S^(f)  represents  the  1-D  spectral  content  of  the 
fluctuating  portion  of  this  property  expressed  in  terms  of  a  unitless  temporal 
frequency  variable  /,  B  is  an  empirically  fit  constant,  and  the  parameter  /o  is 
related  to  the  peak  frequency,  fm.  In  their  original  paper,  Kaimal  et  al.  [1972] 
found  B  =  0.16.  More  recently,  Kaimal  and  Finnegan  [1994]  report  B  =  0.164. 
This  function  is  termed  universal  because  Kaimal  et  al.  viewed  it  as  applicable 
to  any  scalar  property.  Hereafter,  I  shall  refer  to  eq.  (22)  simply  as  the  Kaimal 
curve. 

In  order  to  show  the  appropriateness  of  members  of  the  template  family  in 
approximating  the  Kaimal  curve,  the  Ym  function  1-D  versions  were  fitted  to 
the  asymptotic  behaviors  of  the  Kaimal  curve  at  both  high  and  low  spatial 
frequency.  The  method  of  fitting  involved  shifting  the  transformed  template 
function  vertically  through  multiplication  by  a  constant  and  horizontal  shifting 
by  multiplying  the  frequency  variable  by  a  scaling  multiplier. 

As  an  example,  Ti,  when  converted  to  the  Kaimal  form,  is  expressed  as, 


*i(/o 


0.6/'  (//2  +  6/1 1) 
(1  +  /'2)11/6 


(23) 


This  form  is  fitted  to  the  asymptotic  behavior  using  1.02045  \I/i  (0.47909  /).  The 
resulting  fitted  curves  are  plotted  against  the  Kaimal  curve  in  figure  1  for  cases 
m  =  0  to  3.  As  seen,  the  Ym  functions  exhibit  peaks  which  are  progressively 
narrower  than  the  Kaimal  curve. 

In  order  to  approximate  the  Kaimal  curve  using  these  template  functions,  it  was 
realized  that  two  of  these  functions  could  be  combined  in  such  a  way  that  the 
asymptotic  behavior  remained  the  same,  but  the  peak  region  could  be  better 
approximated. 

Due  to  its  incorrect  limiting  behavior  at  low  frequency  (it  approaches  a  constant 
rather  than  zero)  the  To  function  was  not  used.  Having  ruled  out  use  of  the  To 
function,  the  next  two  members  of  the  T  family  were  used: 


oc  A  Tx  (k,  Co)  +  (1  -  A)  T2  (k,  Cb) ,  (24) 
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Figure  1.  Fits  of  Tm  functions  to  the  asymptotic  behavior  of  the  proposed  curve 
of  Kaimal  et  al.  [1972]. 


Figure  2.  Comparison  between  Kansas  experiment  temperature  fluctuation 
spectra  data  and  Kaimal  et  al.  [1972]  empirical  curve. 


Figure  2  shows  the  data  used  in  this  analysis  along  with  an  empirical  curve  fit 
to  this  data  as  contained  in  figure  20  of  Kaimal  et  al.  [1972]. 

Mathematically,  we  can  relate  this  comparison  to  a  new  form  for  the  outer  scale 


influenced  3-D  refractive  index  spectrum, 


$«(*)  =  PnCl 


>1!  (k£0)2  £a1/3 

r  9i 17/6 

1  +  ( nCa ) 


+ 


(1  -  >li)  (/c£b) 
1  +  (ft  £&) 


4  ~ll/3 

H 

i  23/6 


(25) 


The  £a  and  £b  parameters  appear  as  placeholders  which  are  related  to  the  peak 
frequency  ( fm )  of  Kaimal  et  al.  [1972]  and  have  a  direct  relationship  to  the  outer 
scale,  but  the  exact  nature  of  this  relationship  will  only  be  obtained  as  a  result 
of  the  analysis. 

Using  this  form  for  the  refractive  index  spectrum  we  may  write  the  variance  of 
refraction  index  fluctuations  as, 


OO 

a2n=  I  $n(«)  4tt  k2  dK  =  ^  Pn  C2n  Aa  Ax  (f2Jz  -  ^£2/3  J  +  ^£; 


2/3 

b 


where  A  a  is  an  integration  constant  (see  appendix), 

tt3/233/2 


Aa  = 


nr(|)  r(l) 


(26) 


(27) 


The  1-D  spectrum  for  the  refractive  index  can  be  generated  from  the  3-D 
spectrum  through  another  integration  (see  appendix). 


OO 

Fn  ( Hi )  =2  J  <bn  2tt  Kr  dnr 

o 

[(«!  Fa)2  +  n 


A^-f pnCl C3J3 


+  (1  -A1)^pnC2nCl/3^ 


,11/6 


1  +  (ki  £a)2 

(«i  Fb)4  +  if  («i  Fb)2  + 


(28) 


1  +  (/ci  FbY 
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Spherical  symmetry  of  the  spectrum,  <f>n,  is  assumed  such  that  the  integral 
results  at  negative  Ki  can  be  combined  with  positive  ki  results,  producing  a 
one-sided  spectrum  containing  only  nonnegative  values  of  k\.  Hence,  the  factor 
2  appears  in  front  of  the  integral  sign. 

We  now  combine  the  results  of  the  previous  calculations  to  produce  the 
expression  Tn(/ti)  =  kx  Fn(/ci)/cr2.  Due  to  the  normalization  by  <r2  quantities 
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Pn  and  C\  divide  out.  The  remaining  expression  is  independent  of  the  scalar 
variable  being  studied. 


where, 


tfn(Kl)  £a/3  (KljCa)2+n 


«1 


Ax 


1  + 


(Ki£a)2 


11/6 


+ 


Ai  Ax  r  -  v9i 17/6 


(29) 


1  + 


(«i  Cbf 


Ax  =  A  A 


£“/3  + 17  (ir  - !)  £' 


2/3 

b 


(30) 


This  expression  is  written  in  terms  of  the  variable  /sx,  while  the  Kaimal  equation 
is  written  in  terms  of  the  variable  /.  «x  is  a  spatial  frequency  representation 
along  the  x  axis  which  is  normally  understood  as  the  direction  of  mean  wind 
flow  in  meteorological  studies.  The  /  variable  is  a  dimensionless  temporal 
frequency  variable  related  to  tower  sensed  meteorological  data.  The  two  are 
related  through  the  equation,  K\  z  =  2i r  /,  where  z  is  the  height  above  the 
surface.  To  permit  comparison  let  us  define  the  variable  Ko  analogous  to  /o, 
such  that  kq  =  2n  fo/z.  Eq.  (22)  may  therefore  be  rewritten  in  terms  of  the 
dimensionless  frequency  ratio  ki/k0  =  ///0, 

*"(*i)  B[kq 

ki  1  +  B  (ki/ko)5/3 

The  behavior  of  the  T  («i)  /kx  function  must  thus  obey  B/k0  dependence  at 
low  frequency  and  n1  dependence  at  high  frequency. 

Considering  the  high  frequency  behavior,  cancellation  results  in  an  expression 
for  kq  as  a  function  of  Ai  and  Ax- 


2/3  1 

K°  ~AxAx' 

Rewriting  this  expression,  we  have  the  condition: 

—  =  Ai  (Can o)2/3  +  —  (1  —  Ai)  (£b«0)2/3  • 
At  the  low  frequency  limit  we  obtain  the  expression, 


B_ 

Kq 


n  Ai  £a/3  +  (1  -  Ax)  £ 


72 


,5/3 


Ax  Ax 


(32) 


(33) 


(34) 
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The  denominator  on  the  right  is  just  «0  2^3.  Eq.  (34)  then  becomes, 

jB  =  At  (£„  k0)5/3  +  ^  (1  -  /ti)  (A  «o)5/3  •  (35) 

We  thus  have  two  equations  in  four  variables:  A\,  Ca,  Cb,  and  kq.  However, 
because  of  the  relationships  between  £a,  £*,,  and  k  in  eq.  (25),  and  the  definition 
given  for  the  outer  scale  in  eq.  (4),  we  presume  that  £0,  £5  must  be  proportional 
to  the  outer  scale  length  L0.  k0  appears  to  perform  a  similar  role  in  KaimaPs 
equation,  except  that  it  must  be  inversely  proportional  to  L0.  Hence,  the 
products  £a  «o  and  £(,  k0  may  be  viewed  as  dimensionless  functions  which 
effectively  displace  their  respective  T  functions  along  the  k  axis.  We  then 
have  only  three  unknowns  (7a  =  £a  5  76  =  £b  «o>  and  Ai)  in  two  equations. 
Conceptually,  we  might  view  A\  as  a  free  parameter  which  then  defines  70  and 
76- 

But  there  is  a  complication  that  arises  when  solving  for  an  optimal  value  of  A\. 
From  figure  1  it  is  observed  that  both  Yi  and  T2  lie  above  the  Kaimal  curve. 
This  means  that  the  value  of  A\  may  not  lie  between  0  and  1  since  what  we  wish 
to  accomplish  is  to  subtract  a  multiple  of  the  upper  curve  (T2)  from  a  multiple 
of  the  lower  curve  (Ti),  effectively  pulling  the  Ti  line  downward  in  the  center 
toward  the  Kaimal  curve  and  the  Kansas  data.  This  is  only  possible  if  A\  is 
greater  than  unity.  Unfortunately,  when  A\  >  1  there  is  no  general  solution  to 
both  of  these  equations  simultaneously  for  a  general  B  value. 

The  best  that  can  be  hoped  for  is  to  try  to  approximately  match  the  slope 
behavior  given  in  eqs.  (34)  and  (35),  while  choosing  a  value  for  A\  that  provides 
a  good  fit  to  either  the  available  data  or  the  Kaimal  curve.  Originally,  I  had 
conceived  that  it  would  be  best  to  fit  to  the  Kaimal  curve,  but  this  did  not 
seem  reasonable,  particularly  in  view  of  the  Kaimal  curve’s  performance  at  low 
frequency.  There,  Kaimal  et  al.  [1972]  recommended  a  value  for  B  of  0.16.  Later, 
Kaimal  and  Finnegan  [1994]  amended  this  constant  to  B  —  0.164.  However, 
examining  the  data  it  appears  that  an  asymptotic  constant  of  B  =  0.20  is  a 
better  fit. 

In  any  case,  the  dispersion  of  the  data  is  significantly  larger  at  low  frequencies 
than  it  is  at  high  frequencies.  It  therefore  makes  sense  that  if  only  one  curve 
can  be  solved  for  it  should  be  the  high  frequency  equation.  Therefore,  solving 
for  7b  in  terms  of  Ai  and  7 a,  we  have, 


76  = 


17 

15 


a-M) 


n3/2 


(36) 


For  the  low  frequency  curve  we  have, 


76 


17 
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-1  3/5 


( 1-A 1) 


(37) 
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Of  interest  is  that  we  may  choose  a  value  for  A\  and  then  vary  7a  and  look  at 
a  plot  of  the  function  7b(Al5  7a)  —  i'b{Ai,  7 a)  to  see  how  this  function  varies. 
This  difference  approaches  zero  for  a  narrow  range  of  7 tt  values  regardless  of 
the  value  of  A\.  Thus,  although  there  is  no  solution  to  eqs.  (34)  and  (35)  in 
general,  we  may  approach  close  to  a  solution  and  use  this  result.  A  relatively 
good  choice  of  parameters  that  appears  to  come  close  to  both  the  centroid  of 
the  data  and  fit  for  both  high  and  low  frequency  limiting  behaviors  is  A\  =  8.2, 
la  =  0-44,  75  =  0.5254.  The  7b  result  compares  well  with  Yb  =  0.5221,  though 
it  is  not  exactly  equal.  Figure  3  shows  the  resulting  curve  comparison  with  the 
Kansas  data.  As  can  be  seen,  the  low  frequency  behavior  of  the  proposed  curve 
appears  to  fit  to  the  centroid  of  the  Kansas  data  better  than  the  Kaimal  curve 
of  figure  2. 


Figure  3.  Comparison  between  the  fitting  method  curve  involving  Tm  functions, 
the  Kaimal  empirical  curve,  and  Kansas  data  reported  by  Kaimal  et  al.  [1972], 


With  these  results,  we  may  return  to  the  description  of  the  proposed  3-D 
refractive  index  spectrum  in  eq.  (25).  This  equation  can  be  rewritten  in  terms 
of  fin,  Cl  K0,  7tt,  and  jb: 


M*)  =  ft.C^“/3 


Ai  («7„/ko)  7“/ 
17/6 

1  +  («7a/«o)  j 


(1  -  Ai)  (/C7b//c0) 

7tU/3 

1  +  (/C7b/ /cq)2 

23/6 

(38) 

With  Ai,  7a,  and  7b  fixed,  (3n  a  constant,  and  constant  for  any  given  moment 
and  location,  the  remaining  parameters  are  /c  and  /c0.  Using  the  definition  of 


the  outer  scale  in  eq.  (4) ,  we  find  for  the  parameter  values  listed  above, 

—  « 4.71393;  k0  «  (39) 

Kq  L/0 

Writing  the  same  result  in  terms  of  Ca  and  Cb, 

K0Ca  =  0.4400  — >  Ca  =  2.0741  L0;  (40) 

K0Cb  =  0.5254  — >•  Cb  =  2.4767  L0.  (41) 

The  resulting  outer  scale  spectrum  seen  in  figure  3  accurately  models  the 
refractive  index  spectrum  to  spectral  frequencies  well  below  the  knee  associated 
with  «0  (by  approximately  a  factor  of  5),  thus  significantly  extending  the 
simulation  of  the  refractive  index  spectrum  to  very  low  wavenumbers.  A  plot  of 
the  outer  scale  region  of  the  resulting  3-D  spectrum  is  seen  in  figure  4. 


Figure  4.  Outer  scale  portion  of  the  proposed  spectrum  compared  to  the 
Kolmogorov  and  von  Karman  spectrums. 


In  this  figure  the  newly  developed  spectrum  is  compared  with  both  the 
Kolmogorov  and  von  Karman  spectrums.  While  the  von  Karman  spectrum  rolls 
off  to  a  constant  value  at  low  frequencies,  the  Kolmogorov  spectrum  increases 
to  a  singularity  at  zero  frequency  while  the  new  spectrum  rolls  over  toward  a 
zero  value  at  the  origin.  Of  interest  is  a  comparison  between  the  von  Karman 
spectrum  and  the  proposed  spectrum.  The  von  Karman  spectrum  reaches  a 
maximum  at  a  frequency  slightly  below  the  outer  scale  frequency  k0  while  the 
proposed  spectrum  continues  to  increase  to  a  peak  approximately  twice  that  of 
the  von  Karman  spectrum. 
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2.4  Outer  Scale  Length 


We  are  now  in  a  position  to  discuss  the  association  between  the  outer  scale  length 
and  the  frequency  maximum  measured  during  the  1968  Kansas  experiment. 
In  the  previous  section,  we  connected  the  outer  scale  length,  L0,  to  the  k0 
parameter.  In  this  section,  kq  is  related  to  frequency  maximum  data  collected 
during  the  1968  Kansas  experiment.  This  connection  leads  to  a  definition  of 
outer  scale  as  a  function  of  similarity  parameters  in  the  atmospheric  surface 
layer. 

To  build  this  connection,  note  that  ko  is  defined  as  kq  —  2ir  fo/z.  But  /0 
is  related  to  the  frequency  maximum,  /m,  at  which  Kaimal’s  universal  scalar 
curve  peaks.  This  curve  is  given  by, 


F(f)= 

[l  +  BU7/o)5/3]' 


(42) 


The  peak  of  this  curve  is  found  by  setting  its  derivative  to  zero.  The  location  of 
the  peak  is  denoted  by  fm.  Solving  for  this  peak  location  we  find, 


fo  —  fm 


(43) 


For  the  original  B  value  (0.160),  fo  =  0.2611  fm.  For  the  revised  value,  B  = 
0.164,  fo  —  0.2650  fm.  The  two  are  relatively  close,  and  since  the  Kaimal  et 
al.  [1972]  analysis  was  based  on  the  0.160  value,  we  will  use  this  result. 

We  have  thus  logically  connected  the  outer  scale  value  L0  to  the  frequency 
maximum  location  fm.  We  find  a  functional  form  for  fm  through  an  analysis  of 
Kaimal  et  aVs  figure  18.  The  consolidated  data  from  this  figure  are  included  in 
both  figures  5  and  6  in  this  document.  In  these  figures,  the  abscissa  coordinate 
is  the  ratio  z/L,  representing  the  ratio  of  height  above  the  surface  (z)  to  the 
similarity  parameter  L  known  as  the  Monin-Obukhov  length.  The  first  of  these 
figures  (figure  5)  presents  fm  in  a  linear  plot  format  as  a  function  of  z/L.  The 
second  figure  (figure  6)  presents  fm  in  log  form  as  a  function  of  z/L.  Along 
with  the  data  are  two  curves,  one  for  L  <  0  and  one  for  L  >  0.  The  distinction 
between  these  conditions  is  termed  the  stability: 


fn 


0.460  y+f  +  0.0014  §>0; 

0.036  +  0.0014  f  <  0. 


(44) 
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Negative  L  values  correspond  to  daytime  (unstable)  conditions,  where  parcels  of 
air  displaced  vertically  tend  to  accelerate  away  from  their  level  of  origin.  Positive 


Figure  5.  Comparison  between  the  eq.  (44)  curve  and  frequency  maximum  data 
reported  by  Kaimal  et  al.  [1972]. 


Figure  6.  Comparison  between  the  eq.  (44)  curve  and  frequency  maximum  data 
reported  by  Kaimal  et  al.  [1972]  with  ordinate  axis  plotted  on  a  log  scale. 


L  values  indicate  stable  conditions  where  air  parcels  tend  to  return  to  their  level 
of  origin. 

Although  Kaimal  et  al.  [1972]  included  a  hand  drawn  curve  fit  to  the  data,  the 


25 


above  relationship  is  also  a  reasonable  fit.  This  functional  relationship  reveals 
several  interesting  features  regarding  outer  scale  dependencies  within  the  surface 
layer  atmosphere.  First,  it  reveals  that  there  are  two  pertinent  length  scales  that 
affect  the  outer  scale  length.  Focussing  on  the  results  for  B  =  0.164,  we  find, 
for  nonnegligible  heights,  z, 


(  0.281  V+n:  L  >  0,  \z/L\  >  0.1; 
\  3.592  v/ziT  L  <  0,  \z/L\  >  0.1. 


(45) 


This  functional  relationship  not  only  reveals  that  the  outer  scale  depends  on  both 
the  height  and  the  Monin-Obukhov  length,  it  also  reveals  that  the  dependence 
is  equal  between  these  two  influencing  variables.  However,  the  outer  scale 
exhibits  considerably  longer  lengths  under  unstable  atmospheric  conditions  than 
under  stable  conditions.  Though  this  difference  is  likely  due  to  bouyancy 
differences  between  the  unstable  daytime  and  the  stable  nocturnal  surface  layer 
atmospheres,  the  exact  mechanism  responsible  for  these  differences  cannot  be 
discerned  from  this  study. 

Writing  the  full  form  of  the  fm  function, 


f  0.281  z/  ( yf+^jL  +  0.0030)  z/L>  0; 

Lo  =  <  ,  >  _ _  <  (46) 

[  3.592 zf \y/-z/L  +  0.0389 J  z/L  <0. 

This  form  reveals  explicitly  how  L0  is  dependent  on  the  Monin-Obukhov  length. 
Both  C2  [cf.,  Rachele  and  Tunick,  1994]  and  the  inner  scale  length,  t0  [cf. , 
Tofsted  and  Auvermann,  1991]  are  also  known  to  vary  with  atmospheric  stability 
conditions.  Hence,  modeling  atmospheric  surface  layer  conditions  will  produce 
the  parameters  needed  to  predict  the  form  of  the  refractive  index  spectrum  as 
a  function  of  height. 

2.5  The  Composite  Spectrum 

An  application  involving  this  new  spectrum  is  as  part  of  a  calculation  of 
the  turbulent  coherence  diameter  (chapter  3).  Outer  scale  influences  on  the 
coherence  diameter  can  be  significant  under  certain  circumstances.  However,  to 
evaluate  the  coherence  diameter  it  is  imperative  that  inner  scale  effects  also  be 
included.  Therefore,  it  behooves  one  to  modify  the  refractive  index  spectrum  to 
include  both  inner  and  outer  scale  effects.  Andrews  and  Phillips  [1998]  present 
several  inner  scale  simulation  functions,  including  their  modified  spectrum,  the 
Tatarskii  [1971]  spectrum,  and  the  Hill  spectrum  [Hill  and  Clifford,  1978]. 

These  spectra  will  not  be  considered  here  in  detail.  Rather,  I  note  that  they 
typically  either  ignore  outer  scale  effects  altogether  or  simply  replace  the  k-11/3 
Kolmogorov  dependence  with  a  term  associated  with  the  von  Karman  spectrum, 

(/c2  +  £-2)  ' .  Inner  scale  effects  are  either  represented  as  a  series  of  terms  as 
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in  the  Andrews  and  Phillips  modified  spectrum  or  as  a  superimposed  multiplying 
factor  riding  on  the  main  Kolmogorov  dependence. 

In  either  case,  it  is  difficult  to  analytically  evaluate  the  effects  of  both  inner 
and  outer  scales  because  once  the  spectrum  is  represented  as  a  single  functional 
form  it  becomes  difficult  to  integrate.  Hence,  I  introduce  here  the  composite 
spectrum.  This  new  spectrum  is  composed  of  three  parts.  The  first  two  of  these 
parts  represent  a  summation  of  two  spectra,  one  representing  only  inner  scale 
effects,  the  other  only  outer  scale  effects.  The  high  frequency  portion  of  the 
outer  scale  spectrum  and  the  low  frequency  portion  of  the  inner  scale  spectrum 
both  reflect  Kolmogorov  behavior.  The  summation  of  these  two  spectra  must 
therefore  be  balanced  by  the  subtraction  of  a  Kolmogorov  spectrum  from  the 
sum.  The  spectrum  thus  produced  includes  both  inner  and  outer  scale  effects 
while  simultaneously  permitting  separate  integration  of  the  individual  parts, 
avoiding,  as  much  as  possible,  computational  complexity. 

This  spectrum  is  written  as: 


$n(K,  (.q,  Lq)  —  ^no(^)  Lo)  T  & o )  (47) 


where  is  again  the  Kolmogorov  spectrum  (eq.  (3)),  the  outer  scale 

modified  spectrum,  <f>no,  is  synonymous  with  eq.  (25),  and 


$n/(«) 


0nC^K-n/3FF 


(48) 


is  an  approximation  [Tofsted,  1991]  to  Frehlich’s  [2000]  inner  scale  modified 
spectrum  based  on  an  analysis  of  measurement  data,  where, 


V(x)  =  log10(x)  +  0.72, 


(49) 


Ff(V)  = 


1, 

1  +  0.424  exp  (-2.725  V2)  , 

0.712  [1+ cos  (ft)]  [l  -  0.786  V2] 
0 


V  <  -3; 

—3  <  V  <  0; 
0  <  V  <  0.9; 
0.9  <  V  . 


(50) 


Although  the  value  7.4  appears  in  eq.  (48)  above  as  a  constant,  this  quantity 
is  actually  somewhat  variable.  The  ratio  £0/7.4  is  supposed  to  represent  the 
Kolmogorov  microscale  rj,  but  the  ratio  between  iQ  and  rj  depends  on  the  Prandtl 
number,  which  is  itself  temperature  dependent  [Hill,  1997]. 

An  example  of  the  behavior  of  the  combined  effects  of  inner  and  outer  scale  are 
seen  in  figure  7.  In  this  figure,  the  spectrum  has  been  normalized  by  dividing 
by  the  Kolmogorov  spectrum.  The  inertial  subrange  portion  of  the  spectrum 
thus  evaluates  to  a  constant  value  of  unity  in  the  curve’s  central  section.  The 
inner  scale  bump  section  is  shown  on  the  right  hand  portion  of  the  curve  and 
the  outer  scale  section  is  seen  on  the  left. 
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Though  an  outer  scale  value  of  5  m  was  used,  this  value  is  not  necessarily  typical. 
Given  the  variability  of  outer  scale  on  height  and  the  Monin-Obukhov  length, 
it  is  difficult  to  speak  in  generalities,  yet,  as  is  predicted  by  eq.  (46),  we  should 
generally  expect  the  width  of  the  inertial  subrange  to  be  approximately  an  order 
of  magnitude  wider  under  daytime  conditions  than  at  night  due  to  the  functional 
difference  in  eq.  (45)  between  stable  and  unstable  conditions.  Regarding  the 
inner  scale,  figure  7  highlights  the  positions  of  two  spatial  frequencies  related 
to  the  inner  scale,  £0,  and  the  Kolmogorov  microscale,  77.  The  two  are  related 
through  the  equation,  77  =  4/7.4  [Hill  and  Clifford,  1978],  which  also  explains 
the  appearance  of  the  factor  7.4  in  the  argument  to  the  V  function  of  inner  scale 
adjustment  factor  Fp. 


Figure  7.  An  example  of  normalized  influences  of  inner  and  outer  scale  on  a 
near  surface  refractive  index  spectrum. 


As  seen  in  the  figure,  the  effects  of  the  outer  scale  extend  for  approximately 
one-and-a-half  decades.  Inner  scale  effects  have  a  similar  extent.  Hence,  even 
though  the  inner-to-outer  scale  ratio  is  approximately  a  factor  of  1000,  there  is 
still  only  about  a  decade  and  a  half  associated  with  the  inertial  subrange  for  this 
particular  choice  of  inner  and  outer  scale  values.  The  effects  of  inner  and  outer 
scale  are  thus  on  the  order  of  significance  of  the  inertial  subrange  in  this  case. 

The  question  may  be  asked  whether  this  situation  is  typical.  Studies  of 
atmospheric  effects  on  propagation  often  ignore  inner  and/or  outer  scale  regions. 
The  rationale  for  avoiding  consideration  of  these  effects  most  often  concerns  ease 
of  computation.  In  the  case  of  the  proposed  outer  scale  spectrum  we  have  a  form 
which  is  not  necessarily  easy  to  use,  but  which  is  nevertheless  not  intractable. 
The  main  question  is  whether  such  calculations  are  necessary  and  the  answer  is 


likely  that  under  most  conditions  the  answer  is  no,  for  the  simple  reason  that 
the  outer  scale  may  be  so  large  (because  of  atmospheric  conditions  or  the  fact 
that  the  line  of  sight  is  sufficiently  removed  above  the  surface)  that  other  factors, 
such  as  the  weighting  function  used  in  performing  integrations  with  respect  to 
frequency,  make  the  resulting  calculations  insensitive  to  spectral  dependence  on 
outer  scale.  Nevertheless,  there  will  be  certain  conditions  involving  surface-to- 
surface  propagation  where  both  inner  and  outer  scale  influences  on  the  optical 
turbulence  parameterizations  will  be  significant. 
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3.  Turbulent  Coherence  Diameter 


As  an  application,  consider  evaluation  of  the  coherence  diameter  using  the 
composite  spectrum.  The  coherence  diameter  is  given  implicitly  in  the  equation 
[Hill  and  Clifford,  1981], 

S  oo 

1=4  n2k2  J  ds  J  [i  —  Jo  *n(K,L0(s),CZ(s),i0(s))KdK,  (51) 

0  0 

where  La,  t0,  and  C2  can  all  be  functions  of  path  position.  In  this  equation,  k 
is  the  radiation  wavenumber  ( 2-ir/X ),  s  is  a  path  position  variable,  S  is  the  total 
path  length,  and  r0  is  the  turbulent  coherence  diameter,  often  called  the  Fried 
coherence  diameter  [Fried,  1966]. 

Path  position  dependent  parameters  are  most  generally  considered  height 
dependent.  Since  the  terrain  underlying  the  path  between  object  and  observer 
is  likely  to  be  nonflat  and/or  a  slant  path,  the  values  of  these  quantities  will  vary 
in  nonlinear  fashion.  However,  at  a  given  moment  for  a  given  point  along  the 
optical  path,  the  governing  parameters  can  be  considered  constant.  It  is  thus 
possible  to  directly  evaluate  the  inner  integral. 

To  evaluate  this  equation  we  begin  by  introducing  the  new  length  variable, 
X  =  sr0/{2.\S).  Using  the  composite  spectrum  introduced  in  the  previous 
chapter,  we  may  write  eq.  (51)  as, 

Po 

1  =  4n2  k2  —  [ pn  C2(XS/Po )  X5/3  [G(X,  Lo)  +  H(X,  l0)  -  ft]  dX,  (52) 
Po  J 
0 

where  the  path  integration  variable  has  been  replaced  by  X,  pa  =  rQ/ 2.1, 
and  G,  H,  and  unity  represent  separate  results  for  different  components  of 
the  turbulence  spectrum.  G  represents  the  combined  outer  scale  effects,  H 
represents  the  combined  inner  scale  effects,  subtraction  by  the  constant  ft 
represents  the  Kolmogorov  effect,  and  La  and  lQ  may  vary  as  a  function  of 
range. 

To  solve  for  ra  involves  an  iterative  method  which  requires  an  initial  value  of  ra  to 
begin  the  search  procedure.  In  general,  rQ  will  be  greater  than  this  initial  value, 
which  is  normally  computed  based  on  a  Kolmogorov  spectrum  assumption. 
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3.1  Kolmogorov  Spectrum  Results 

For  the  Kolmogorov  portion  of  the  spectrum,  the  inner  integral  may  be  written, 

OO 

(3nC2  j[\- J0  (kX)]k~11/3  KdK  *  0.033  x  1.11833  C2  Xb/3.  (53) 

o 

The  constant  terms  then  combine  as,  47r2  x  0.033  x  1.11833  «  1.45695.  Through 
further  analysis,  we  find  that  the  integration  constant,  1.11833,  is  actually  the 
approximation  of  the  expression, 

**  =  25/3  r(  ].  1  /6)  r  (1 1/6)  ”  1  11833440°-  («) 

Returning  to  the  main  double  integral  in  eq.  (51),  we  have, 

s 

1  =  4jr2  k2  Pn  N  j  Cl  (|T| )5/3  ds,  (55) 

0 

which  can  be  solved  explicitly  for  ra  as, 


2.1 


-\  -3/5 


(47r2/?nN)  k2  S  J  C2(Sw)w5/3dw 


(56) 


where  w  =  s/S  is  a  dimensionless  path  position  variable  and  C2  may  vary  as 
a  function  of  range.  This  form  corresponds  to  typical  integral  representations 
of  the  coherence  diameter  for  the  Kolmogorov  refractive  index  spectrum  [cf., 
Roddier,  1981]. 

One  important  factor  to  consider  in  the  use  of  this  equation  regards  the 
associated  path  position  weighting  function,  w5/3.  In  imaging  problems,  this 
weighting  assumes  the  imager  is  located  at  range  S  (w  =  1).  In  propagation 
problems  (laser  beams),  this  weighting  function  is  normally  reversed  such  that 
the  laser  source  is  located  at  the  origin.  In  this  case,  a  weighting  function  of 
(1  —  id)5/3  is  used. 

3.2  Outer  Scale  Spectrum  Results 
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To  evaluate  the  outer  scale  cases  we  now  introduce  three  new  dimensionless 
variables,  Ya  =  Ca/X,  Yb  =  £b/X,  and  u  =  kX.  Given  the  spectral  form 


associated  with  the  outer  scale,  evaluation  of  the  inner  integral  in  eq.  (51)  will 
involve  the  solving  of  equations  having  the  form, 


OO 

Gm(c/x)  =  x~5/3 /[i- MkX)}  rm(K,  c)KdK 


oo 

=  Gm(Y)  =  Y11/3  J  [1-Mu)} 


[i  +  (My)2]m+11/6 


udu. 


Thus,  the  factor  X5/3  sifts  out  of  these  equations  just  as  it  does  from  the 
Kolmogorov  expression. 

After  considerable  math,  it  can  be  shown  that, 


OO 

G i  (Ya)  =  yo11/3  J  [l-Jo(u)] 


— (^ya)2 17/6  u  du  =  ya5/3 

[i  +  («y0)2]17/6  55 


is  21/6  r  i 


55  r  (I)  Ya/6  L2 


0^1/6  Ur  -yaK5/6 


It  can  also  be  shown  that  in  the  limit  of  large  Ya, 


OO 

Gi  (Ya)  =  Y^3J  [1-  Jo(tt)] 


(uYgf 
[1  +  (uYa)*f 


.dutt#  —  —Y  x^3.  (59) 

55 


For  the  second  function  we  find  the  following  results: 


oo 

c«(n)  =  n11/3/  [i  -  Jo(«>] 


216  21/6  f  23 


W)4  udu  =  ™y°/° 
(o,VA 2i23/6  935  6 


[1  +  («y6)2]23/6 


935  r  (|)  Yb1/6  L24 


In  the  limit  of  large  Y, i,  for  this  function  we  find, 


OO 

G2(Yb)  =  Ybll/3  J  [1-JoW] 


(uYb)4 

[i  +  (u  n)2]23/6 


udu^x-  —y: 


In  these  results,  the  functions  Kv(x)  are  modified  Bessel  functions  of  the  second 
kind,  order  v  [cf.,  Kreyszig,  1972].  However,  the  raw  integration  results  are 
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expressed  in  terms  of  modified  Bessel  functions  of  the  first  kind,  I„(x),  which 
are  related  to  functions  of  the  second  kind  through  the  relation, 

Ku^  =  2  sin(mr)  ~  '  (62) 

These  calculations  also  involve  use  of  the  recurrence  relationship, 

2  jj 

Iu(x)  =  Iu-2{x) - Iu-i(x), 

X 

and  the  series  expansion  form  of  Iu(x), 

™  x2m+v 

^  ^  Z_/  22m+u  m\  F(m  +  v  +  1)  ’ 

Combining  these  effects,  we  generate  the  function, 

G(Ya,  Yh)  =  A1  Gi(Ya)  +  (1  -  A,)  G2(Yb). 

Let  Y  =  La/X .  At  large  Y  we  have, 

G(Y)  «  N  -  0.78078  Y~1/3. 

The  general  solution  is  plotted  in  figure  8.  In  this  figure  the  approximation 
function  was  used  for  Y  values  greater  than  100  due  to  numerical  stability 
problems  encountered  in  evaluating  the  Bessel  functions  at  high  argument  values. 
The  approximation  at  this  level  is  accurate  to  5  digits  of  accuracy. 


(63) 

(64) 

(65) 

(66) 


Figure  8.  Outer  scale  influence  function  for  the  coherence  diameter,  plotted  with 
respect  to  the  parameter  Y  =  L0/X. 


3.3  Inner  Scale  Spectrum  Results 


The  results  for  the  inner  scale  effects  must  generally  be  evaluated  numerically 
due  to  the  form  used  in  describing  the  inner  scale  modifier  to  the  spectrum.  But 
the  form  of  this  calculation  is  essentially  similar  to  the  one  made  for  the  outer 
scale, 


OO 

H(£0/X)  =  X" 5/3  J[  1  -  J0{kX)]  k"8/3Ff 


V(^ 

\  1A, 


dn 


(67) 


This  calculation  can  be  written  as  a  function  of  the  dimensionless  variable  Z  = 
iD/X.  The  results  of  these  calculations  are  shown  in  figure  9. 


Figure  9.  Inner  scale  influence  function  for  the  coherence  diameter,  plotted  with 
respect  to  the  parameter  Z  =  iQ/X. 


To  study  the  behavior  of  this  function  at  large  Z  we  perform  the  change  of 
variables  u  =  kX, 


H(Z)  = 


11.4 /Z 

f  [l-JoM]u-8/3f> 


0 


du. 


(68) 


The  upper  limit  on  this  integration  is  due  to  the  upper  cutoff  of  the  bump 
spectrum  function  Fp  combined  with  the  definition  of  the  V  function.  For  large 
Z  the  limits  of  the  integration  will  become  small,  to  the  point  that  the  quantity 
[1  —  Jo(u)]  u-8/3  can  be  approximated  by  w-2/3/ 4,  using  only  the  first  two  terms 
of  the  Bessel  function  expansion. 
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Considering  the  behavior  of  this  portion  of  the  integrand,  and  considering  that 
the  Fp  function  is  almost  constant  throughout  its  nonzero  range,  it  should 
perhaps  come  as  no  surprise  that  the  integral  evaluates  approximately  as  3u1/'3/4 
with  an  upper  limit  proportional  to  1  /Z.  Analysis  of  calculated  integration 
results  for  101  <  Z  <  103  yields  a  functional  form  of  approximately, 

H(Z)  »  0.61068  Z~1/3,  Z>1.  (69) 

This  result  is  used  to  evaluate  the  H (Z)  at  Z  values  where  numerical  techniques 
fail  due  to  integration  difficulties  near  the  singularity  at  u  =  0. 


4.  Conclusions 


In  this  report,  the  modeling  of  outer  scale  effects  on  optical  turbulence  within 
the  atmospheric  surface  layer  has  been  considered.  We  first  discussed  the  means 
used  in  treating  outer  scale  effects  by  previous  investigators.  We  considered 
prior  limitations  in  the  definition  of  outer  scale  and  in  the  use  of  the  standard 
von  Karman  spectrum.  A  new  definition  of  outer  scale  was  given  such  that 
a  particular  point  in  the  energy-containing-eddy  region  of  the  refractive  index 
spectrum  is  used  to  denote  the  outer  scale  wavenumber,  ft©.  A  development  was 
also  described  which  relates  outer  scale  effects  on  the  refractive  index  spectrum 
to  the  research  results  of  Kaimal  et  al.  [1972]  and  a  subsequent  reanalysis 
by  Kaimal  and  Finnegan  [1994].  These  results  link  the  outer  scale  value 
to  frequency  maximum  data  for  temperature  fluctuations  in  the  surface  layer 
atmosphere.  Analysis  of  this  data  and  the  subsequently  proposed  dimensionless 
height  ( z/L )  parameterization  of  the  frequency  maximum  (/m)  indicates  that 
the  outer  scale  is  not  merely  a  function  of  height  above  the  surface,  as  has 
been  commonly  assumed,  but  rather  is  a  function  of  height  and  the  Monin- 
Obukhov  length,  firmly  establishing  that  the  outer  scale  is  a  function  of  both 
of  the  length  parameters  associated  with  the  surface  layer.  As  part  of  this 
analysis,  a  new  outer  scale  influenced  spectrum  was  developed  which  avoids 
the  zero  limit  problem  of  the  von  Karman  spectrum.  Lastly,  the  implications 
of  the  outer  scale  were  considered  in  combination  with  inner  scale  effects 
within  a  proposed  composite  refractive  index  spectrum.  The  resulting  spectrum 
permits  the  separable  evaluation  of  inner  and  outer  scale  effects  on  atmospheric 
propagation.  As  an  application,  we  considered  the  evaluation  of  the  turbulent 
coherence  diameter  based  on  this  new  spectral  model. 
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Appendix 

Integration  Methods 


In  section  2.3.3  of  this  report,  several  integrations  were  performed.  In  this 
appendix  the  methods  used  in  performing  these  integrations  are  described.  The 
first  of  these  operations  is  related  to  the  Ym  family  of  functions  introduced  in 
section  2.3.2.  In  2.3.3,  these  functions  must  be  integrated  into  one-dimensional 
(1-D)  form.  Figure  1  contains  a  family  of  such  integrated  curves.  Second,  we 
integrate  these  functions  to  generate  the  variance  statistic  cr£. 

These  integrations  can  all  be  performed  using  the  3.259.3  result  from  Gradshteyn 
and  Ryzhik  [1980]  (hereafter,  G&R): 


00 

J  xx_1  (1  +  axp)~ M  (1  +  / 3xp)~ v  dx  =  ^  oTxtp  B 
o 


x  2Fi  (  v,  fi  +  w,  1  -  - ]  ; 
VP  a/ 


(A-l) 


where  |  arg  a|  <  7r,  |  arg  /3\  <  7r,  p  >  0,  and  0  <  5RA  <  2 +  v).  The  function 
B(x,  y )  is  the  beta  function, 


B(x,y )  = 


r(*)  r(y) 


T(rr  +  y)  ' 

The  2-P1  function  is  Gauss’  hypergeometric  function,  given  as, 

,  o  \  ,  oc(3  a  (a  +  1)  P  ((3  +  1)  2  , 

JF1(a,/3i7;2)  =  l  +  —  z+  ^^jTjyTTTJ  2  +' 


(A-  2) 


(A  —  3) 


In  performing  the  evaluation  of  <s\  begin  by  considering  the  integral, 

OO 

a2  =  J  4>n  (/«)  47 r  n2  dn.  {A  -  4) 

o 

Removing  constant  terms  from  the  integral,  we  have  two  resulting  integrations 
to  perform: 

/  H  +  7 aJ7 dK = &  <^-5) 

J  [1  +  (KCa)2]  J  [1+U2] 


where  u  =  nCa,  and, 

[1  +  (k£j,)2] 


7  ("A)4£;i/3  ,.2d,._rv3  J 
J  [1  +  ucmwk  dK~Cb  J 


vr 


[1  +  U2 


2i23/6 


du , 


(A -6) 
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where  u  —  nCt,- 

We  can  write  these  two  integrals  in  the  more  general  form, 


OO 

/(»)  =  / 


(i + x2),,+r,/6 


(A -7) 


This  integral  can  be  solved  using  the  cited  G&R  formula  where  the  following 
substitutions  are  introduced:  a  =  0=l,  p  =  2,  X  =  2n+1,  fx  =  u,  p  +  u  = 
n  +  5/6.  All  the  components  here  are  real,  and  the  final  condition  is  met  where 
2n  +  1  <  2n  +  5/3.  Hence,  we  obtain  the  result, 


IM  =  1  r(»+l)  rfl) 

()  2  r(n+|)  ' 


(A  -  8) 


Here,  the  hypergeometric  function’s  main  argument  (1  —  P/a)  evaluates  to  zero, 
and  hence  the  function  itself  evaluates  to  unity  and  drops  out  of  the  final 
expression. 

To  evaluate  the  transformation  from  three-dimensional  (3-D)  spectra  to  1-D,  we 
begin  by  considering  the  integral: 


OO 

2  J  +  «2^  27 TKrd,Kr. 


(A -9) 


Again,  removing  constants,  we  find  the  following  two  integrals  at  the  heart  of 
the  problem: 


[1  +  *?  C2a  +  K2rJCl]17/6 


ACy’  dt\>tp 


OO 

=  C ^  f 


(u2  +  tx2) 

[1  + 1>2  +  u2]17/6 


idu,  (A  —  10) 


where  u  =  KrCa  and  v  —  KXCa\  and, 


'  (4  +  4?  cl  c"* 
[1  +  Kf  c\  +  kI  cif/e 


Kr  d,Kr  =  C^3 


(v2  +  U 2)2 

[1  +  v2  +  u2]23/6 


udu\  (A  —  11) 


where  u  —  KrCb  and  v  —  nxCb-  To  simplify  these  integrals  we  can  remove  the 
v2  dependence  from  the  denominators: 


(v2  +  u2) 

[l  +  u2  +  «2]17/6 


udu  = 


[1  +  v2] 


.2117/6 


[1  +  A 


,2117/6 


(A  - 12) 


[1  +  v2]17^6  J  [1  +Au2] 


2117/6 
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where  A  =  (1  +  v 2)  1;  and, 


(v2  +  u2)2  _  v4  r  u 

[1  +  u2  +  u2]23/6  W  “  “  [1  +  u2]23/6  J  [l  +  ittf3/6 


[1  +  u2]23/67  [1+^m2]23/6 


du  {A  -  13) 


OO 


[1  +  ^2]23/6  7  [l  +  ^U2j23/6 


where  ^4  =  (1  +  u2)  *. 

These  forms  are  thus  seen  to  devolve  to  the  solution  of  the  following  general 
integral: 


(l  +  /lx2)n+5/6 


(A  -  14) 


where  0  <  m  <  n  —  1.  This  form  can  also  be  solved  using  the  previously 
mentioned  G&R  result.  Here,  we  use  A  =  2m  +  2,  p  =  2,  a  =  /3  =  A,  p  =  v, 
and  p  +  v  =  n  +  5/6.  In  general,  we  find  a  result, 


_ _ ix  -  1  d  ,  „2)(-+l)  +  1)  H(n  -  m  -  1/6) 

(i +^xt+5/6  2  *  ^  u  n„+5/6)  -  (j4  15) 
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